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Abstract. The troublesome multi-center molecular integrals of Slater type orbitals (STO) 
in molecular physics calculations can be evaluated by using the Fourier transform and proper 
coupling of the two center exchange integrals. A numerical integration procedure is then 
readily rendered to the final expression in which the integrand consists of well known special 
functions of arguments containing the geometrical arrangement of the nuclear centers and the 
exponents of the atomic orbitals. A practical procedure has been devised for the calculation of 
a general multi-center molecular integral coupling arbitrary Slater-type orbitals. Symmetry 
relations and asymptotic conditions are discussed. Explicit expressions of three-center one- 
electron nuclear attraction integrals and four-center two-electron Coulomb repulsion integrals 
for STO of principal quantum number n=2 are listed. A few numerical results are given for the 


purpose of comparison. 



I. Introduction 


Because the calculation of three and four-center integrals is one of the most complicated 

problems involved in the quantum theory of molecules and solids' , it has been a stumbling block 

to molecular ab initio calculation. The nature of the complications is due to the complicated 

angular dependence associated with different centers and the presence of square roots of radial 

dependence in describing relations between functions of different distances. In the past, two 

general procedures have been employed. The method of Barnett and Coulson* is based on a one- 

center expansion in which the numerical work commonly is very tedious and prone to error. The 

other method is based on the adoption of Gaussian orbitals either via expansion’ or 

transformation' . Although the Gaussian type orbitals are computationally convenient, they are 

nevertheless inadequate because of poor radial dependence, e.g., the nonfulfillment of the cusp 

condition as well as the rapid falloff at long distances. Therefore, numerious Gaussian 

orbitals must be introduced in order to produce an accurate result. In the past, a few 

systematic investigations of molecular two center integrals over Slater-type orbitals (STO) led 

to general and explicit, but formal, solutions , 5 - which commonly cannot be used efficiently 

for numerical calculations. Recently, instead of tackling the STO, which supposedly give a more 
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accurate electronic distribution, Filter and Steinborn introduced the so called "B functions" 
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and "A functions"” which are closely related to the reduced Bessel functions and Laguerre 
functions, respectively. Linear transformations exist between the STO and these two types of 
functions. A general expression' ' has been given for the three-center nuclear attraction 
integral and the multicenter repulsion integral over the "B functions". However, the advantage 
of these functions is dubious both because the expression itself is still complicated, involving 
multiple numerical integration, and because additional effort is required to perform the 
summation to obtain the STO integral values. We have already suggested 1 2 that, by using the 
Fourier transform technique and properly coupling a pair of two-center exchange integrals, the 
multicenter molecular integrals can be handled in a straightforward manner. For example, a 
four-center two-electron Coulomb integral can be cast into a simple triple numerical integration 
with finite domain. As the complicated geometrical arrangement is connected vectorially and 
appears in the argument of the functions in the integrand, a simple numerical procedure can be 
directly applied. However, we recently noted that the polynomials appearing in the integrand 
are related to the modified Bessel functions. 1 5 Therefore, by means of their recurrence 
relations, the Final compact expressions are readily obtained. 


II. Basic integrals and procedure 


A set of exchange integrals which are commonly applied to the atomic collisions of heavy 
particles has the following form 


ii 2 i 2 

I (a,b,c,d) = f r" 1 r” 2 (P, -r )(p 2 .r )..(p r ) 

n, ,1, ,n, ,1, J A B A A 1, A 


(q, .r B )(q, r B )-(q, - r B ) ex P(- c r A - d r B +ia r A +ib r B ) dr 


( 1 ) 


where, as before' 2 p, ,p 2 ..p and q, ,q 2 ..q are arbitrary constant vectors and, in 

*1 4 


particular, can be set to unit normal vectors along the Cartesian axes; a and b are real 

constant vectors. The constants c and d arc real and positive. Integers n, and n, * I and 

1 and L > 0. The vectors r , r are the radius vectors of the electron measured from nucleus A 
1 ^ A B 

1 4 

and B respectively. It has been shown that Eq. (1) can be expressed as a simple one- 
dimensional integral such as 


I 

n, 


I. +1, 


A 


(a,b,c,d) = (-i) 


(P, 


•V-Ci 


•V' 


.v b )..«,| 


•V 


n, -1 n 2 -1 

(-3/3c) (-3/3d) I (a,b,c,d) 

1 >0 ,1 »0 


( 2 ) 
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where I (a,b,c,d) -2n f 

1 *0 9l 10 J . 


, exp(i (J. R - A.R) 


du 


X = u(l-u)la + bl + uc +(l-u)d 


R = r -r 

A B 


P = ua-(l-u)b, 

Also, the following differential operations can be easily established 


(3) 

(4) 

(5) 


V X = V A, = u(l-u)(a+b)A 

a b 

a X/dc = ucA , i) X/dd = ( 1 -u) d/x 


V (p.R) = uR , V (p.R) = -(l-u)R 

a b 

d(P.R)/a c =o, a(p.B)/ad=o 


( 6 ) 

(7) 

( 8 ) 
(9) 


In adopting this expression, we suggested' 2 that the multi-center molecular integrals of 
STO can be constructed from these two-center exchange integrals. Although in the final 
expressions numerical integration still has to be performed, the troublesome multi-center nature 
is connected in a vectorial relation which appears only in the argument of the integrand. Thus, 
after the geometry and the orbitals are given, the integrals can be evaluated in a straight- 
forward manner. The method has been found to be expedient, and accurate results have been 
obtained. In this communication, we note that the parametric differentiation in Eq. (2) is 
commutative and associative. Therefore, if the solid harmonic operators are introduced, in 
relation to Eq.(l ), we can further define another set of two-center exchange integrals such as 

J . , (a,b,c,d)= f r ■ ' Y, (rJr*Y, (rj 

n, ,1, ,m, ;n 2 ,1 3 ,m 2 J A 1 , m, A B 1 2 n\ B 
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n . -1. -1 n, -1, -1 , , . .. . , 

r A r B exp( - cr A dr B +,ar A + 'br B )dr 


( 10 ) 


lnEq.(10)n, ,1, ,m, and n, ,1 2 ,m, are the quantum numbers of the orbitals. In connection 


with Eq. (2), we obtain 


1, +1, 


J . . (a,b,c,d) =(-i) 

n, ,1, ,it>, ;n 2 ,1 2 

n, -1, -1 n 2 -1 2 -1 

(-3/3c) (-3/3d) I 


\ m , ( VWV 


(a,b,c,d) 


2 50 52 50 


*. H * 


n. -J, -1 


n 2 -4 -1 


=(-i) y, (V ) y (V ) (-3/3c) 
1 , m, a 1, m 2 b 


(-3/3d) 


2ncdf 1 e lp R R k 2 (A,R)u(1-u)A du 
() 


(ID 


where k is the modified spherical Bessel function. Furthermore, the solid harmonic 
n 

1 6 

gradient operator y (V ) introduced can be interpreted as in the Condon-Shortley 
1, >m, a 

1/2 

convention , for example, y = -(3/8 ir) (3/3 a +i3/3a ) etc. For practical purposes, the 

i *i x y 

differential operations have to be performed under the integral sign in Eq.(l 1) which is 

straightforward but laborious. However, with the aid of Eq.(4)-Eq.(9) and the recurrence 

relations of k , the procedure can be greatly simplified. The technique of solving the 
n 

multicenter molecular integral by coupling a pair of two-center exchange integrals has already 
been suggested. To illustrate the idea and to demonstrate the simple operations, we write down 


the three-center nuclear attraction integral, denoted as L, and the four-center two-electron 
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Coulomb repulsion integral, denoted as M, coupling some STO, the numerical solutions of which 
are known to fib difficult.. ■ 


We Write 


L = fy (£ A ,r )- V , ^n’ r n^ dr 

J n, ,1, ,m, A A r n, ,l 2 ,m, B B 


( 12 ) 


Using the identity and referring to Fig. 1, we have 

1/r = 1/2 it Je^A^AcVq dq 

Eq.(12) can be written as 

L = C (t) C (£) f e K| R AC 12k q dq 
n. An, B J 


1, +1, 


n, -1, -1 


n j -4 - 1 


(■') K - a/J V , ( - 3/3 V WWro'V 

2 "«A«B/;- ei|1 '* AB R AB K - (XR AB )/X ' U(1 - U)dl "a= q ,b=0 


where p = ua -(l-u)b, X =u(l-u)la+bl +u£; +(l-u)t 

A B 

After performing the parametric differentiation, the angular parts of q are integrated 


(13) 

(14) 


analytically by utilizing the plane wave expansion exp[iq.(uR -R^ )], and the radial parts of q 
and u are integrated numerically. The four-center two-electron Coulomb integral is defined as 


follows: 

M 


m !\ ,1, *,«A* A>\ ,1, .m, ~ \ .1, , m , V 1 C> 


12 
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( 15 ) 


y . (£ n ,r, )dr,dr, 

n 4 ,1 4 ,m 4 D D 


In Eq.(l5), for example, the indices n, ,1, ,m, are the quantum numbers of the orbital; r, 
denotes the radius vector of electron number 1 referred to an atomic nucleus B as origin; ^ is 
the exponent parameter of STO etc. 


We replace 1/r, , by its Fourier transform 

1/r, , =1/2* J e iqT ' 2 /q dq =1/2* J e' q(R AC +ri C' r ' A ) /q dq (16) 


Notice that expressing r, , as R +r, -r, is quite arbitrary and offers no advantage 

AC C A 

over other choices as long as they form a vectorial relation. 


M = C (%. )C ( % )C ($ )C (4 )/(2* q ) 
n, A n, B n, C n 4 U 


1 , 


* A 


1, 


1 , 


* A 


1 , 


/ r ’ A Y |, mf r ' A )f| B Y |. .n <r 'B )r, C Y l, m, (r ’C )r, D \ J. r ' D } 

n, -l, -1 n.-4-l n,-l, -1 .j iqr iq.R,„ 

r ' A r ' B r 'C r ’D e Ae Ce AC 

'^A r ' A'^B r ' B ^C r ’ C'^D r ’ D 


dr, dr, dq 


(17) 


By using Eq.(10), we write Eq.(17) as 

M = C (£)C (L)C (^ r )C (t ) J e ,q ' R AC /(2rc q ) 
n, A n, B n, C n 4 D J 


dq 


1 , + 1 , + 1 } + 1 4 


n, -1. -1 


n, -1, - 1 


(-0 


K-OT5 a ) <-wy y,, m , <v a )y l, m, (V b ) 


^ ^ R AB (XR AB )A ’ U(1 - U) ' a =-q,b=O dU 


8 



n, -1, -1 


n 4 'I, - 1 


[ (-d/d$ c ) 


(-d/^ D ) 


y l, m , ( V y l. V 


( 18 ) 


i 



iP’.R 


« CDR CD K - <rR AB )A ' U ' <IV)l a'= q ,b’=O d "' 


where P and A. are as defined before and 

p’=uV-(l-u’)b\ X* =u’(l-u’)la , +b’l I +u’$* +(l-u’)5* D (19) 

Because the angles are generated via parametric differentiation with respect to the 
constant vectors a, b, a’ and b\ it is necessary to invoke the conditions, a = -q, b = 0, a’ = 
q, b’ = 0 at the end of the differential operations. Then, the integrand consists of a finite 

I 7 

sum of angular orbitals in variable q. By properly coupling the orbitals and using the 
plane-wave expansion exp|iq.(R -uR +u’R )], we can integrate analytically the angular part of 

AL A15 LU 

q. Then the radial parts of q , u and u’ integrations are calculated numerically in a 
straightforward manner, such as by Gaussian quadrature of other numerical means. 

111. Symmetry relations 


By inspection, there are simple symmetry relations that can be readily employed; in this way the 
number of independent integral expressions derived is reduced to the minimum. 
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Three-center Nuclear Attraction Integrals 


Referring to Fig. 1, suppose we have an expression (see Appendix A ) for the following 
integral. For brevity, we write i= n, ,1, ,m, , j= rij ,1, .m, etc. 

If the expression of the following integral is known 

V I vw !yw r 


( 20 ) 


we consider the following integral 

L i.r/V 5 A’ r A ) f c 7i (5 B' r B ,dr 

where L\ .denotes the integral for which the nuclear center C has undergone an 
J.i 

inversion operation, i.e., if r =lr - R I, then r = lr + R I, Furthermore, through the 


inversion of the electronic coordinate which is equivalent to the operation taking r to -r^, 

r to -r 4 , Eq. (21) becomes 
B A 


( 21 ) 


L j.r/V^A'- r B ) f c V i (? B-- r A )dr 
=( - ,) ' /VVB ) L V i ( ^B' r A )dr 


1, 


=(-l) 


1, 

+1: 

=(-l) 



i, j V ^A o 

The interpretation of Eq.(22) is that the value of integral L . can be obtained by 


( 22 ) 


taking the complex conjugate of the value derived by means of the expression of L.’ . in which 

bj 
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the exponents are properly interchanged with the inverted nuclear coordinate C. Finally, a 
phase factor is multiplied to obtain the integral value. 

Four-center Two-electron Coulomb Repulsion Integrals 


Referring to Fig. 2, suppose we have 

'ij;k,f K ( V' A^/V' B> C ,V I ( V- D )dr ' dr - 

and are interested in evaluating another integral of the following form 

'k,l;i,j = / VV' A> V 1 ( V' r, , VV’ C^/V’ D )dr ' dr - 


The integral value of I . is, in general, different from that of I. . Because 
the subscripts 1 and 2 are dummy indices, Eq.(24) can be easily rewritten as 

Vl;i,f / V^A’ r ' A> W’ B } r , , VV' C> V j ( V' D )dr ' dr “ 

■W AMC - BMD - 4 A" t C 5 B" t D 1 

In Eq.(25), exchange of the nuclear center is implied, e.g., A to C, C to A etc., 

including the exchange of exponents in the expression of I. . Essentially, the nuclear 

1>J jK,l 

center is reordered. 


(23) 


(24) 


(25) 
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Next, consider the following integral 


£ 1 jjc 

Vyjt- / y j < V' A> V i ( V' B> t , , V i ( V’ C> \ ( -D' r - D> dr ' dr > 


Comparing to Eq.(23), we immediately find that 


( 26 ) 


I. .. =1. .. ,(A — B, C < — * D, ^ ♦ » £ , £ « * ^ ) 
j,i;l,k i ,j;k,l A B C D 


(27) 


IV. Asymptotic Behavior and Nuclear Center Collapse 


AB 


With this formulation, more information of the asymptotic behavior of the integral is 

easier to obtain, at least in a qualitative sense. For example, it is obvious that the 

integrand decreases exponentially as R^ and increase linearly through k^. However, if R 

and R are kept finite, the integrand decreases in a dampened oscillatory fashion with respect 

to R through j . In another words, the integral value varies to a smaller extent with respect 

to a change of R than of R and R , as we would expect. The formula is equally valid if the 
AC AB CL) 

nuclear centers are collapsed into three, two or one center by simply taking the nuclear 

distance to zero, e.g., R = 0, or R =0, (three center), or R =R =0 (two center), or 

Ad CD AB CD 

R =R =R =0 (one center), respectively. We notice that k (x) goes to infinity at the origin, 
Ad LU AL n 

but of course, this does not necessarily mean that the integral is not finite. However, for n > 


2, there is a much stronger singular point at the origin; the straightforward procedure by 
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directly taking e.g., R -* 0 would certainly lead to some numerical problems. Fortunately, by 


m 


close inspection, we find that, at all times, the term R k (&R), m > n+1 appears in the integrand. 


Using the series representation of k 


n 

-x, M (n+1)! , -1 (n+2)! . -2 

K n x = e /X| ' + I! r(„) <2x) + 2T r(n-l) 2x + 1 

-x " 1 -k 

= e /x Z (n+ - , k) (2x) 

0 ^ 


where the Hankel’s symbol (n, k) = - 


F( 1/ 2+n+k) 


k! T( 1/2+n-k) 


(28) 


we therefore, obtain 


1 i m R m k (XR) =0 , m > n+1 
R->() " 


L2n >! _ 

,.n.n+l . 
2 X n! 


, m =n+l 


(29) 

(30) 


i.e., for the collapsed center case, Eqs.(29) and (30) can be used in the integrand to avoid any 
numerical difficulties. 


V. Numerical Example 


In order to test this method, a simple Gaussian quadrature integration scheme has been applied 
to the derived formulas for the two electron multicenter integrals for the methane molecule. 

The results are listed in Table I along with the well known published results of Shavitt and 
Karplus 4 , and the results computed using the Barnett-Coulson 2 method. We used the identical 
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parameters and the same notation as reported by Shavitt and Karplus, which are summarized in 


Table I. 


Table I. Some multi-center-two electron integrals for methane 



X 

y z 

Orbital 

Symbol 

Exp. 

c 

0.0 

0.0 0.0 

ls(C) 

c 

5.700 

H, 

0.0 

0.0 -2.0 

2p(C) 

x,y,z 

1.625 

H, 

1.8856180 

0.0 0.6666667 

Is(H.) 

h. 

i 

1.000 

H, 

-0.94280904 

1.63299316 0.6666667 




H< 

-0.94280904 

-1.63299316 0.6666667 






Results, according to number of quadrature points 




This method 

Shavitt-Karplus 

Barnclt-Coulson 


Integral 

12x12x16 

20x20x24 12x12x16 

20x20x24 



(h, h, ;h, h, ) 

0.03068148 

0.03068213 0.03068196 

0.03068187 

0.03068247 


(h, h, ;h, h, ) 

0.03569343 

0.03569397 0.03569398 

0.03569368 



(h, h, ;h 2 h, ) 

0.09569741 

0.09570521 0.09570384 

0.09570556 



(c h, ;h, h, ) 

0.01298606 

0.01274316 0.01286213 

0.01271293 

0.01274093 


(x h, ;h, h, ) 

-0.00752578 

-0.00752560 -0.00752573 

-0.00752559 

-0.00752559 


(z h, ;h, h, ) 

-0.04765937 

-0.04766222 -0.04766145 

-0.04766277 

-0.04766263 
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(c I), ;c h, ) 
(c h, ;x li, ) 
(c h, ;z h, ) 
(x h, ;x h, ] 
(x h, ;/. h, ) 
(y h, ;y h, ) 
(z h, ;x h, ) 
(z h, ;z h, ) 
(c c ;h, h, ) 
(x c ;h, h, ) 
(z. c ;li, h, ) 
(x x ;h, h, ) 
(x 7 ;h, h, ) 
(y y ;h, h 3 ) 
(z z ;h, h, ) 
(x x ;z h, ) 
(x z ;x h, ) 
( 7 . z ;/, h, ) 


0.01191481 

0.01121649 

0.01191464 

0.01121651 


0.02017340 

0.01981324 

0.02018282 

0.01980944 

0.01980903 

0.00602142 

0.00588692 

0.00602480 

0.00588562 

0.00588522 

0.03557018 

0.03556931 

0.03556692 

0.03556880 

0.03556858 

0.00 149034 

0.00148956 

0.00148892 

0.00148937 

0.00148940 

0.02231921 

0.02231975 

0.02231903 

0.02231976 

0.02231959 

- 0.07790808 

- 0.07791435 

- 0.07791088 

- 0.07791394 

- 0.07791452 

0.00855243 

0.00854907 

0.00854768 

0.00854864 

0.00854849 

0.16644702 

0.16654179 

0.16864892 

0.16656439 


0.00236735 

0.00237118 

0.00236900 

0.00237002 

0.00237004 

- 0.00167397 

- 0.00167668 

- 0.00167514 

- 0.00167586 

- 0.00167581 

0.14192876 

0.14192778 

0.14192799 

0.14192784 

0.141927 

- 0.00143553 

- 0.00143592 

- 0.00143598 

- 0.00143598 

- 0.00143579 

0.13676011 

0.13675834 

0.13675888 

0.13675864 

0.136758 

0.14091369 

0.14091243 

0.14091270 

0.14091244 

0.140912 

- 0.22143208 

- 0.22142793 

- 0.22142658 



- 0.01062253 

- 0.01062202 

- 0.01062207 



- 0.25503541 

- 0.25502346 

- 0.25502024 
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VI. Conclusion 


Using the integral transform technique, we have derived explicit expressions for one- 
elcctron three-center nuclear attraction molecular integrals coupling the n=2 atomic Slater 
orbitals which are listed in Appendix A. Using Eq. (20) and similar expressions given in 
Appendix A, we can write the expressions for two-electron four-center Coulomb repulsion 
integrals, which are considered among the most difficult molecular integrals to evaluate. It 
is neither necessary nor feasible to list all the expressions; however, we give in Appendix B, a 
representative set. These expressions contain well behaved and well known special functions to 
which numerical procedures can be easily applied. Although the development is parallel to what 
was reported earlier! 7 the recognition of the modified Bessel functions in the integrand and the 
gradient operators introduced allow us to write directly the integral expression coupling any 
arbitrary orbitals in a much more compact form. As pointed out earlier, the merit of this 
approach lies in the fact that the complicated geometrical nature of a multi-center integral is 
properly handled by means of the vectorial relations and is manifested in the integrand through 
the argument of special functions in a natural way. Therefore, it may have better numerical 
behavior than most aforementioned methods. For a small set of molecular integrals, we have 
encountered no numerical difficulties as clearly borne out by our numerical example. However, a 
more thorough investigation will be conducted in the future. By means of a simple symmetry 
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argument, the number of independent integral expressions can be greatly reduced. We also have 
shown the asymptotic behavior of the integral with respect to the internuclear distance and 
expect no numerical difficulties if provisions are properly taken to treat the collapsed center 
case. With perseverance, other integral expressions involving either higher orbitals or other 
operators can be equally well deduced. Despite the fact that this method cannot be easily 
generalized, i.e., explicit expressions have to be developed for each particular integral, this 
method, nevertheless, offers indeed a simple and compact solution to this notoriously difficult 
molecular problem. With further refinement and efficient numerical procedures, we believe this 
approach will make a contribution to the field of molecular calculation. We also point out that 
this formulation is quite suitable for a small computer, as was demonstrated in our numerical 
example. A typical run of a set of integrals on a IBM PC AT is of the order of a few minutes. 
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Appendix A 


Three-center nuclear attraction integrals 

We list only the integrands with the argument of the functions omitted wherever possible, 
he full expression should be interpreted in the following sense, for example: 

L = <2p 0 A - lsB> 

r c 

| 4 2 2 

= (,^ b C, (5 a )C, « b )/7 l-j.(q»)M».R AB )R AB u (1-u) <| 


Y, 0 (co)A + j 0 (q«) k 2 AR ab ) R ar u (1-u) Y, 0 (R An ) A Idu dq 


AIT 


We further substitute q = tan(v) to use v as an integration variable with the limit from 
0 to n/2. In the above equation, we have 

to = i|R ab* r ac ’ r ab =r b' r a’ r ac =r c' r a 

2 2 2 2 

and X is the positive root of X = u(l-u)q + u£ +(l-u)£ B 

C (^) = (24) (2n+1)/2 r (2n)! 
n 

U =<lsA-‘ lsB> = ^ A ^ B C 1 (^ A )C,(^ B )u(l-u)j 0 R An K 3 A A (Al) 

« 4 2 2 3 5 2 

U=<2sA- lsB> = £ Q (^. )C, (E )j. [ R.u (l-u)t k, A -R.-uO-u)*, A |A 
r B A B Ali A All 

(A2) 
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5 2 2 2 2 


.,=<2sA- rc 2sn> = C 1 (5 A )C,(^)j.|R AB „ (l-») s a 5 b . 4 a -R ab «,A 

2 2 3 2 

(u^-Kl-u)^ ) u ( 1 -u)+R ab u(1-u)k 2 A 1/jc 
L« = <2p„ A i IsB > = 2 /vA^ b Q (^ a )C, (^R^u* (l-u)f-j, k, (l-u)q 

Y, 0 («)A +j„ k 2 Y. 0 (R ab )A 1 


L, =<2p ±1 A- 1 lsB> = 2/y7i^ B Q (4 a )C, (^)» (l-u)R AB |-j, k, (l-u)q 


* A 1 


* A 2 


Y, ±1 (o>)A +j 0 k 2 Y, ±) (R ab )A 1 


L. = <2p 0 A J 2sB > = 2/VnC 2 (^)C 2 (^ A [- R AB u 2 (1-u)’ k, q j, ^ Y, . (w)A 

+R AB j° Kj u (1_u) ^n Y| » ( r a J a +u ^- u > r aJ. q 


’B 


AB 


AB 


Y o (co)A -u (1-u)k 2 R ab Y, 0 (R AB )j 0 A ] 


4 =<2p ±1 A 1 2sB> = 2/7jiC 2 (^ a )C 2 ^ b >^ a I- r ab u 0- u ) k '< q j. ^ 0 Y , ± , (w)A 
+ R An j u u (1-u) $ Q Y, ±1 (R ab )A +u (1-u) K, R A J, q 


AB 

Y* (®)A -u (1-u)k 2 R Y + (R. n )j 0 A ] 

i i| Ad j ±i Ad 


AB 


* A 


L« = <2p 0 A - 2p 0 B> = 44 a ^ b Q (^ a )C 2 (^ b )[k 4 q R Ag u (1-u) A (j 2 (1/5 jc) ' 
Y 2 , (co)-j 0 14* )-j, k, R AQ q u (1-u) (2u-l)Y, 0 (R^Y, 0 (w) 

3 <223 

A +3j 0 Kj R u (1-u) A 14 k 


(A3) 


(A4) 


(A5) 


(A6) 


(A7) 
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(A8) 


-jo R ab u (1-u) Y, . (R ab )Y, . (R ab )A I 


L, =<2p 0 A- 1 2p ±1 B>=4^ A ^ B C 2 (^ A )C J (^ B )[q ^ R AB u 0-«) A (3/20rc)' 4 j 
Y 2 ±1 (to)-j, K, R AR qu (1-u) (uY, 0 ( R A g)Y, ±1 (co)- 


(l-u)Y, 0 (to)Y, ±1 (R ab ))A 

-jo Ki R ab u 2 (1-u) 2 Y* 0 (R ab )Y, ±1 (R ab )A | 
1 


5 2 3 3 4 


L, „ =<2p ± , A- 2p ±1 B>=4 ^ a ^ b C 2 (^ a )C 2 (^ b )[k 4 R ab q u (1-u) A 
(-j 0 /4 jc-(1/20jc)' 4 j 2 Y 2 0 (£>)) 


+j, R* AB qu (i-u)* {(i-u)Y, ±1 ( r ab )y, + ± , A) 

-uY, ±1 («)Y* ±I (R ab ))A -j. R ar u 1 (1-u/ tc 2 
Y ,* ± , A ab )Y, ±I (R ab )A 2 +3j 0 k, R ab u (1-u/ A /4k] 


(A9) 


(A 10) 
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Appendix B 


Four-center two-electron repulsion integrals 

We list only the integrands and suppress the argument of the functions wherever possible. The 

full expression should be interpreted in the following sense, for example: 

M, 0 = <2p , A IsB - 2sC lsD> 
i r , - 


- ^B C ’ <V C ' <V C ' <V C ' W'fj. S. '- j ' (<1W> K - aR AB> R 


AB 


* A 


u (1-u) qY Jw)/X +j 0 (qft>) k 2 (XR AR ) R AB u (l-u)Y + ( R An )A 1 


>± 


±, AB 


[R CD < (VR^K O-OW’ R CD< 1 dudu’dq 

We further substitute q = tan(v) to use v as an integration variable with the limits from 0 to 
tc/ 2. In the above equation, we have 

“= r ac- uR ab +ii ' r cd- r ab=W r ac= r c 'V r cd= r d r c 
2 22 22 22 2 
X = u(l-u)q + u^ A + (1 -u)E, b , X' = u’(l-u’)q +u’^ c +(l-u)^ D 

C (^) = (2^) (2n+1)/2 V(2n)! 

n 

M, = <lsA IsB - 1 IsC lsD> =C, (^C, (I^C, (^ C >C, (^ D )^ B ^ D j. /( 2lt ) 


u(i-u)k 2 rx R AB u’(i-u’K R CD /r 


M, 


<2s A 1 sB - 1 sC 1 sD> =C 2 JC, ^ R )C, (t, ^C, D ) j B /(2 k) 


(Bl) 


[ r ab k 3 u (l-u^^A -R ab k 2 u(1-uK b A 1 [u’(1-u’)^ c ^ d < r cd A’ I 


(B2) 
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M, = <2sA IsB - 2 S Cl S D>=C,(5 4 )C,(yC,(^C,(^)j„ /(2k) 


1 2 

2 2 


= D 


J 5 


|r ab k > “ 0-«>6 a V^ -R ab «. "0 -u) 5 b A ] 

( r cd k i' »' (l-«')^c^D /r ' R CD K; 1 

M. = <2sA 2sB ^ ^ IsC lsD> -C, (£, A >C, (ij^C, « D > j. /(2 k) 


2 2 2 4 4 


|r ab k - u {1 - u) 5 a$b a - r ab k ' (u V (N u) V 

u(l-u)/X +R ab u(1-u) kJX ]!u’(1-u’)5 c $ d <R£ D A’’] 

M 5 = <2sA 2s B ^ 2sC lsD> =C 2 (^)C 2 (^)C 2 (^C, (^) j 0 /(2 k) 


2 2 2 4 4 


| r ab k 4 u (1-u) § a ^ b A - r ab k »(u? a +(1-u)^ b )u(1-u)A +r ab u (1-u)k, A 1 

r R CD K » u ’ ( 1 - u ’W r - r cd< i 

M 6 = <2sA 2sB ^ 2sC 2sD> =C 2 (^)C 2 (^ R )C 2 (^C, (^) j 0 /( 2 k ) 

[ R AB K< U °* u) ^A^B A ' R A B Kj (U ^A +(1 ‘ U) ^B )U(1 ' U)A +R ab 11 0-u)k 2 A j 


2 2 2 4 4 


,r cd*>’ (1 - u,) Vd a ’ - r cd< (u V (1V V 

u’d-iOA’ +R cd u’( 1- U ’)< A’ 2 i 


M, = <2p 0 A IsB - ^ IsC lsD> = C 2 (^C, (^C. (^C, 

A 5 


(B3) 


(B4) 


(B5) 


(B6) 


3. 2 


K qR AB j, u (l-u) Y, 0 (C0)A +j 0 tc 2 R ab u (l-u)Y, 0 (R ab )A Mu’d-u’X r cd A’ 1 


(B7) 


M, = <2p^_, A IsB^ ^ IsC 1 sD> = C 2 (^)C, (i^C, (i^C, (^aVcV^ 
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[K, qR AB j, u (1-u) Y, ±1 (w)A +j 0 tc 2 R aq u (l-u)Y, ±1 ( R A0 )A Jfu’(l- 
1 ' ,)k; r cd a ’ 1 

M, = <2p 0 A IsB - 1 2sC lsD> = C 2 (^C, (^ fi )C 2 (l^C, 

[«c, qR AB j, u (1-u) Y, o (w)A +j 0 k 2 R ab u (l-u)Y, 0 ( r ab )A 1 
[R cd k; u’ 2 (1-u’)^ d A’ 5 - R qd < u’(1-u’)^ d A’ ] 


M, 


= <2p ±1 A IsB - 2sC lsD> = C 2 (^C, (l^C, (^ A 


[K, q R AQ j. u (1-u) Y, ±1 (<o)A +j 0 K 2 R ar u (l-u)Y, ±1 ( R A g)A 1 

ir cd k; u ’ 2 (1 - u,) Vd a ’- r cd< 1 

M, , = <2p 0 A IsB - r 2sC 2sD> = C 2 A >C, (^ R )C 2 (^)C 2 
Ik, q R A0 j, u (1-u) Y, 0 (w)A +j u k 2 R ar u (l-u)Y, # ( r ar )A 

|R cd k; u ' (1 ' u ’ )2 44 A>< r cd K j ’ ^ c + (i-u’K d ) 

3 3 2 

u’(l-u’)A’ + r cd u ’(1-u’)k’ A’ I 

M, 2 = <2 P±1 A IsB - 1 2sC 2sD> = C 2 (y)C, (yC 2 (^)C 2 ^ R /^ 

[k, qR AB j, u (1-u/ Y ; ±I («)A +j. k 2 R ab I.' (1-u)Y* ± 1 (R ab )A I 

I r cd k: u ’ (1 - u ’ )2 44 - V - r cd k ’ <“ VoW d > 

3 3 2 

u’(l-u’)A’ +R cd u’(1-u’K A’ 1 

M, , = <2p„ A 2p 0 B i IsC lsD> = (^)C 2 (yC, (yC, (y 


(k 4 r ab q u (1-u) A (j 2 (1/5 jc)' Y 2 0 (co)-j 0 /4 k) 


(B8) 


(B9) 


(BIO) 


(Bll) 


(B12) 
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5 2 2 


+j, K Afi q u (1-U) (uY, 0 (R ab )V, „ (w) 


4 2 2 3 


-(l-u)Y, 0 («)Y, 0 (R ab ))A +3j 0 Kj R aq u (1-u) A /4rt 


5 2 2 


-j. R aq u (1-u) Y, . (R ab ))Y, . (R ab )A 1 u'(l-u’K R cd /X' 

M, , = <2p 0 A 2p ±1 ' IsC 2 S D> = 21; ^C, (^)C, (^>0, (^C, « D ) 


I 2 

S3 3 4 


fq R ab u (1-u) A j 2 (3/2 Ok)' A Y 2 ±1 (©)+ 


5 2 2 


j, R AQ qu (1-u) (uY, . (R ab )Y, ± , ((oHl-u)Y, 0 (<o)Y, ±1 (R ab »A 


5 2 2 


-jo k 2 R ab u (1-u) Y, 0 (R ab )Y, . (R ab )A ] 

[R cd < u ’ (1 - u, )^ c 4 d A’ -R cd < u*(1h, V j 

M, s = <2p ±1 A 2p ±1 B ^ 2sC 2sD> =21^^ (^)C 2 (^)C 2 (^) 

5 2 3 3 4 


K R AQ q u (1-u) A (j® /4 jc+(1/20jc)' A j 2 Y 2 0 (a,)) 


5 2 2 


-j, k, R AQ qu (1-u) {(l-u)Y, ±1 (R AB )Y* ±t (a) 


3 5 2 2 


-uY, ± , (co)Y* ±I (R ab ))A -j, R ab u (1-u) k- 2 


* A 


4 2 2 3 


Yl ±« (R ab )Y| ±' (r ab )A +3ju Kj r ab u (1 ' u) a ,4n] 


2 2 2 4 4 


ir cd k;u- a-,,-) t c t D A- R cd < («■ VO-’Kd) 

3 3 2 

u’(l-u’)A’ +R cd u’(1-u’K A’ ] 

M, . =<2p. A 2p. B J ^ 2p„ C 1 S D>«45 a ? b 5 c ? d C, (^)C, (4 q )C, (l^C, (5^* 


4 3 2 


5 4 3 3 3 


K r cd A ’ u ’ 0- u ’> (-*. r ab a u (1-u) q {(27/28)’ A (1/5 tc) j, Y, 0 (w) 
-(9/20jt)j, Y, 0 (w))+k 5 R ab /X q u (1-u) (2u-l){j 0 /4ic-(l/5ic)' A j 2 Y 2 . (©)} 


(B13) 


(B14) 


(B15) 
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Y, „ (R AB )+q j. Y, 0 Mu (1-u) { k, R ab A Y, 0 (R ab )Y, 0 *(R ab ) - 3k, R^a’ !4k } } 

-< R CD A> u> (1 ' u)Y ' » R AB A u (1 ' u) q /4n ^ Y * » .<“) 

(1/5 re)' 4 } -k, R ab A q u (1-u) (2u-l)j, Y, 0 (R^JY, 0 M+{k 2 R^/X Y, 0 (R ^) 

A 4 J 2 2 


Y ' 0 ( r ab ) - 3,c > r ab a /4ji )u (1 ' u) ^ (B16) 

M, 7 = <2p ± , A 2p ±1 B - 1 ^ 2p ±l C lsD> =4^^^, (^)C 2 (^ fi )C 2 (l^C, (^Ai 

4 32 2 543 33 / 

K R CD A ’ U> (1V) { ' K%R AB A U (1 " u) q ( -(3/20jt)(2/7)' 2 j } 

^ 5 3 2 *2 2 A 

Y, ±1 (w) -3/1 Ore j, Y, ±l M} + k, R^A q u (1-u) { uY, ±1 (R AB ) 

a Yj 0 (w)(1/20tc)’ 4 +j 0 /47t)+(l-u)j 2 (3/20 k) 4 Y* 2 (w)Y, (R . _)} 

± ± AB 

-qj. Y; ±t Mu (1-u) {k 2 R ab A Y, ±1 (R AB )Yr ±I (R AB )-3tc, R ab A /4 k ) 

< R cd A’ u’ (1-u’)Y* ±1 (R cd ){k 4 R ab A u (1-u) q {j 0 /4 tc +j 2 Y 2 0 (©) 

(1/207C)’ 4 } -k, R ab A qu (1-u) (uY* ±I (R^Y, ±1 M+(l-u)Y* ±1 (w)Y, ±1 (R^) 

5 1 * A A 4 5 2 2 


+ {k 2 R ab A Y, ±1 (R ab )Y, ±1 (R ab )-3k, R ab A /4k }j 0 u (1-u) )| (B17) 

M, . = <2p 0 A 2p 0 B - r ^ 2p 0 C 2p 0 D >=^ A ^^ C ^ (S A )C 2 (^)C 2 (^)C 2 (^k 

543 3 54334 / 

K R CD A ’ u ’ O' 11 ’) KR A b A 11 O-u) q {( 3 / 357 c)(lA)' /J Y 40 (u»j 4 
-(9/\4k)(\/5k)' 4 j 2 Y 2 0 M+(9/80 tc )j 0 } +k, R^A u* (1 -u)' 
q Yl 0 ^ AB ^ 2u - 1 ) ^ ( 1 / 5 tc)( 27 / 28 )' 4 j, Y, 0 M-( 9 / 20 jt)j, Y, 0 (w)} + 

q {K > r ab a y ' 0 (R ab )Y| 0 (R ab ) 3,c ’ r ab a /4k} 1,2 (111)2 

{jo /4 tc -(1/57C) 1 4 j 2 Y 2 o (co)} | 


27 



' k ’ r cd a ’ u ’ 2 (1 ' u ) ^ r ab a u (, ‘ u) {(2uM)Y ' o(R AB ){(l/5rt) 

(27/28)' 4 Y, 0 («)j 3 -(9/20jt)Y , 0 («)j, } }+k 3 R^A u (1-u) q Y, 0 (R^) 

Y, o (R CD )(2u-l)(2u’-l)(j 0 /4* -0/5*)' 4 Y 2 . j 0 )+(k 2 R^A Y, 0 (R^) 

Y, 0 (R ab )-3*c, R ab A /4 k)u (1-u)' q(2u’-l)j 1 Y, 0 (®)Y, 0 (R^)] 

+u’ 2 (1-u’/ {< R cd A’ 2 Y, 0 (R cd )Y, . (R cd )-3k; R cd A’ /4k}[k, R^A 
u (1-u) q {j 0 /4rt -(1/5 jc)' 4 j 2 Y 2 0 (w)} - k, R ab A u (1-u) q(2u-l) 

Y, „ (R ab )Y, . (o,)j, +{ k 2 R ab A Y, 0 (R ab )Y, . (R ab >3k, R ab A /4k)]) (B18) 

M, , = <2 P±1 A 2p ±1 B I" 2p ±I C 2p ±l D> =«$ A $ B * C * D C, (^)C 2 (^)C 2 (^)C 2 (^)tc 


$ 4 3 


3 4 1 ) 4 




{<R A’ u’ (1-u’) [k 4 R ab A u (1-u) q { (3/140jc)(l/rc)‘ 2 Y 4 0 (<o)j 

/ A * J * * * 

+(3/14tc)(1/5jc)' 2 j 2 Y 2 0 (co)+3j 0 /(40 jc )} -k, R^A u (1-u) 


* A 


4 3 


q* { (l-u)Y, ±1 (R ab ){-(3/20tc)( 2/7)‘ 4 j, Y* ±1 («)-(3/10n)j, Y,' ±1 («)} 
uY* (R A g) { -(3/20rc)(2/7) i 4 j, Y, ±1 (®H3/10it)j, Y, ±I (w)} + 

q { (1/2071)’ 4 j 2 Y 20 («)+j 0 /4tt}{K 2 R ab A Y* ±1 (^ AB )Y ' ±' (R A b)' 3iC>R AB A /4rt}1 

532 2 54333 /\ 

+ K ; R cd A’ u’ (1-u’) lK 4 R ab A u (1-u) q (-(l-u’)Y, ±1 (R^) 

{(3/20rc)(2/7)' 4 Y* ±1 (w)j 3 +(3/l()7t)Y* ±1 (w)j, }+u’Y* ±l (R^) 

{(3/10tc)Y, x , (w)j, +(3/20 tc)( 2/7)' 4 Y 3 +1 (w)j, ))-{k 2 R^A 


* A 


+ ' 

3 2 


Y *±' ±' ( r ab »-. r ab a /<4,t))u <l u) jl |Y ' i' (R CD )Y *±' ^ 


A * A 


5 3 2 


(l-u’)-Y, ±1 Ay; ±1 (R cd )u’) + k 3 R A q [-(1-uXl-iOY. ±I (R^Y, ±I (R^) 
(3/IOjc)' 4 j 2 Y* ±2 (w)-{u’(1-u)Y, ±1 (% b ) y .* ± . (R cd ) 
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+u(l-u’)Y* ±l (R ab )Y, ±1 ( R CD )}(j. /4n +(1/20 k)‘ U Y 2 0 (w)j 2 ) 
-uu’Y* ±1 (R AB ) y * ± . (R CD )(3/10icy A Y 2 ±2 (»)j 2 ] +{k 2 R cd A 


* A 


Y. ( R ^) Y . ( r ^>3k; /( 4 ")) u ’ 0-“’) I k < R ab A u (1 ‘ u) 


±* CD ± CD 1 CD' 
q 0„ /4 tc +(1/20ic)' A Y 2 0 (w)j 2 )-k, R. r A q u (1-u) j 


* A 


AB 

A * A 


{ ( 1 -u) Y , ±1 (R ab )Y, ±1 (fl))+uY, ±I (<o)Y, ±l (% B ))+jo 

<*■ r ab a Y * ±- (Ar ab )Y - ±- (R ab } - 3 ^ r ab a /(4lt)m 


M, 


1 


0 = <2p ±1 A 2p 0 B 2sC 2sD> = 2C 2 (^)C 2 (^>C 2 (^C, (^^ 


5 4 3 


* A 


{q j 2 (3/207t)' A Y* ±1 (w)k 4 R^A u (1-u) -qj. [Y, ±1 (w)Y, 0 (R AB )0-«) 


* A 


A 2 


2 2 * A 


(B19) 


-uY, ±1 (R ab )Y, . (o»]u (1-u) Kj R ab A -k 2 R ab J 0 u (1-u) A Y, ±1 (R^Y, „ (R AR )} 
[ r> cd < u>J a-“’) 4 44 A * - R CD K; 


u’(l-u’)A’ + R u’O-u’X A’ 1 


(B20) 


M 2 , = <2p ±1 A IsB r 2p ±1 C lsD> =2C 2 (^ A )C, (^Q (^C, (^ A ^ C ^ D 


4 2 2 3 


4 3 2 2 2 / * A 

K’R u* (1-U’) A* {K,R a q U (1-u) (3/1071)' 2 Y 2 ±2 (co)j 2 


2 * A ♦ A 


-qj, k 2 R u (1-u)A Y, ±1 (R ab )Y, ±1 (w))+k 2 ’ R’^A’ u’ (l-u’){qj, k, R ab A 


AB 

2 * A * A 


2 * A * A 


U (l-u) Y, ±l (R^)Y, x , (co)+j 0 k 2 R ab u (1-U)A Y, +l (R An )Y, +1 (R rn ))l ( B21 ) 


±‘ AB 1 +' ' CD 7 


±’ v CD 7 1 ±' 

M 2 2 = <2p 0 A IsB - 1 2p 0 C lsD> =2C 2 (^)C, (^)C 2 (^^C, (^ D )^ B ^ D 

[k’ R u’ (1-u’) A’ (k',R A q u (1-u) (jo Mic-j* (l/5ic)' /l Y 2 0 (a,)) 

CL) AL> > 

-qj, k 2 R at ,u (l-u)A Y, 0 (R aq )Y, 0 (w)}+k 2 ’ R’ nn A’ u’ (l-u’){qj, k, R ar A 


AB 


CD 


AB 
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(B22) 


1 A 


u (1-11) Y, 0 (R cd )Y, 0 («)+j, R ar u (1-U)A Y, . (R ab )Y, 0 (R cd )>I 
M 2 , = <2p ±1 A 2p 0 B - 1 ^ 2p ±1 C lsD> =8 C 2 (^)C 2 (yC 2 (^C. 


4 3 2 2 3 


3 4 3 3 


(-k; r cd A’ U* (IV) [q (1/4^X18/35)* A Y, ± , (i)Y 0 . k 4 R A u (1-u) j, 


5 3 2 2 


/ * A A 

vl / 2 


,/* 


-q h K, R ab A u (1-u) {(6/2071)' 2 Y 2 ± , (<o)Y, 0 (R AB )(l-u)-u(3/207t) 

Y * ± . (Ar ab )Y -±' ( " )} ' qY *±' ( “ )Y » • J* K * r ab a u (1 ’ u) y *±> (R ab )Y ‘ • (R ab : 


4 2 2 

^CD y 


-<R A’ u’ (l-u’)Y* (lL n )Y 00 { -q j 2 k, R A u (1-u) (3/20^)* A 


±' CD *° 0 ' M J2 *'AB' 

Y, ±) (co)+qj, k, R ab A u (1-u) f Y, ±1 («)Y, 0 (R ab )(1-u)-uY, ±1 (R^Y, 0 (©)} 

+J. ^ R ab A u (1-u) Y,* ±1 (R ab )Y, . (R ab )J 

M 2 4 = <2p ± , A 2p 0 B ^ ^ 2p ±1 C 2p 0 D> =8C 2 (^ C 2 (^ fi )C 2 (^ C )C 2 (^k 
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